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ABSTRACT 

Context. One of the main aims of the ESA Rosetta mission is to study the origin of the solar system by exploring comet 
67P/Churyumov-Gerasimenko at close range. 

Aims. In this paper we discuss the origin and evolution of comet 67P/Churyumov-Gerasimenko in relation to that of comets in general 
and in the framework of current solar system formation models. 

Methods. We use data from the OSIRIS scientific cameras as basic constraints. In particular, we discuss the overall bi-lobate shape 
and the presence of key geological features, such as layers and fractures. We also treat the problem of collisional evolution of comet 
nuclei by a particle-in-a-box calculation for an estimate of the probability of survival for 67P/Churyumov-Gerasimenko during the 
early epochs of the solar system. 

Results. We argue that the two lobes of the 67P/Churyumov-Gerasimenko nucleus are derived from two distinct objects that have 
formed a contact binary via a gentle merger. The lobes are separate bodies, though sufficiently similar to have formed in the same 
environment. An estimate of the collisional rate in the primordial, trans-planetary disk shows that most comets of similar size to 
67P/Churyumov-Gerasimenko are likely collisional fragments, although survival of primordial planetesimals cannot be excluded. 
Conclusions. A collisional origin of the contact binary is suggested, and the low bulk density of the aggregate and abundance of 
volatile species show that a very gentle merger must have occurred. We thus consider two main scenarios: the primordial accretion 
of planetesimals, and the re-accretion of fragments after an energetic impact onto a larger parent body. We point to the primordial 
signatures exhibited by 67P/Churyumov-Gerasimenko and other comet nuclei as critical tests of the collisional evolution. 
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1. Introduction 

Promoting the quest for the origin of the solar system is funda¬ 
mental to the concept of the ESA/Rosetta mission. A comet nu¬ 
cleus - specifically, that of 67P/Churyumov-Gerasimenko (here¬ 
inafter 67P) - was chosen as the target because of the common 
notion that comet nuclei are icy planetesimals formed beyond the 
snow line in the nascent solar system and the evidence that their 
atomic and molecular compositions provide a good match to the 
primordial material out of which the solar system was built. In 
particular, comet 67P was selected as the back-up target after the 
failure to launch Rosetta toward 46P/Wirtanen, partly because it 
was expected to be relatively fresh after having recently encoun¬ 
tered Jupiter (in 1959) and thus had its perihelion distance much 
reduced. 

The formation of the icy planetesimals is currently viewed 
within either of two different scenarios. One is the classical 
model of hierarchical accretion (Weidenschilling l20081l . and the 
other considers the assembly of gravitationally unstable “pebble 
clouds” by gas-grain instabilities in the solar nebula (Johansen et 
al. 120141 1. Recent simulations (Wahlberg Jansson and Johansen 
12014b have shown how planetesimals in the size range from 10 
to 1000 km can thus be produced. 

Independent of which scenario is relevant to comet 67P, it 
seems clear that its ultimate origin has to be sought in the same, 
very early environment. The exact place is more difficult to spec¬ 
ify. The region inside the initial orbits of the giant planets can 
likely be discarded, since cometary ices would not condense 
there. Because of the gas drag presented by the solar nebula, km- 
sized planetesimals formed in the giant planet zone are unlikely 
to be gravitationally scattered outward in order to populate the 
well-known cometary reservoirs (Brasser et al. 12007b . Thus, the 
trans-planetary space, hosting a massive disk of planetesimals, 
is currently favored as the birth place of comet nuclei. 

The issue of collisional evolution among those planetesimals 
has been discussed for almost 20 years. Davis and Farinella 
(11997b modeled the collisional evolution of the Edgeworth- 
Kuiper Belt (EKB) and argued that its members with radii of 
a few km - like typical Jupiter family comets - are multigen- 
erational fragments formed by the splitting of larger objects. 
Stern and Weissman (12001b considered the formation of the 
Oort cloud out of planetesimals formed in the giant planet zone. 
They showed that, during the course of this gravitational scat¬ 
tering process, comet nuclei would be destroyed by collisions 
with small debris. However, following the introduction of the 
Nice Model for the evolution of the solar system (Tsiganis et al. 
12005b . the EKB is no longer thought to be the direct survivor 
of the initial comet populatiorQ, and according to Brasser and 
Morbidelli (12013b . the actual Oort cloud formation may have lit¬ 
tle to do with the classical concept of scattering of planetesimals 
from the formation zone of the giant planets. 

In this paper we discuss some of the evidence that OSIRIS 
imaging of the 67P nucleus may bring regarding its origin and 
evolution. We argue that the bi-lobate shape is likely due to 
a gentle merger of the two lobes as separate bodies in a low- 
velocity collision. This conclusion is supported by additional 
observations, such as the presence of thick layering. We place 
this into the framework of the collisional evolution issue as it 
currently stands, and we highlight the importance of resolving 
the remaining issues about the fundamental question, whether 
the 67P nucleus is the result of an early merging of two distinct 


* The cold EKB may indeed be a remnant of the primordial disk, but 
it carries a very small fraction of the total mass (Fraser et al. 12014b . 


planetesimals or was re-accreted in connection with a major col¬ 
lision long after the planetesimals were formed. Meanwhile, we 
also point out that, in any case, this nucleus does not appear to 
have completely lost its primordial signatures and may thus re¬ 
tain its significance as a witness of planetesimal formation in the 
nascent solar system. 

In Sect. 13 we discuss the role of past erosion in shaping the 
67P nucleus and evidence regarding its bi-lobate shape and the 
interpretation as a contact binary. We summarize the evidence 
concerning structural features of the nucleus and the origin of 
the contact binary in a low-velocity collision. Sect. [3 gives an 
estimate of the collision rate in the trans-planetary disk, which 
illustrates that it is generally unlikely for km-sized planetesimals 
to survive intact during the expected lifetime of this disk, and in 
Sect. |4] we discuss possible ways for comets to escape destruc¬ 
tion in the event that 67P should prove to be such a survivor. Our 
conclusions are summarized in Sect.|5] 

2. Overall shape and structural features of 67P 

2. 1. Previous erosion 

Thermo-physical models reveal that 67P could have lost a sur¬ 
face layer of up to several hundred meters thickness due to the 
accumulated activity (Sierks et al. I2015I Keller et al. l2015b dur¬ 
ing its previous orbits in the inner solar system. This estimate 
is based on the median number of such orbits as estimated from 
long-term integrations for hctitious objects similar to 67P, and 
the specific evolutions of the perihelion distance of the comet, 
which largely governs the erosion rate. Figure [1] shows an exam¬ 
ple of the large variety of previous evolutions of the perihelion 
distance of 67P that may lead into orbits close to the present 
one according to long-term integrations. In general, these evolu¬ 
tions have brought the comet’s perihelion closer to the Sun, but 
Figure [1] covers only part of the whole time spent by the comet 
in the inner solar system. In many cases, the comet had made 
earlier visits into orbits with small perihelion distance, which 
also contributed to the cumulative erosion. We note that the dy¬ 
namical time scale found by these integrations may be consis¬ 
tent with the capture time scale recently estimated by Guzzo and 
Lega (12015b . although further work is needed to confirm this. 

While the use of the current shape and spin axis orientation 
cannot have a drastic influence on those results, the extrapola¬ 
tion of the current activity level may possibly be an important 
source of error. In particular, there is a concern about the pos¬ 
sibility that the comet has undergone long periods of dormancy 
(Kresak ri987b . when its gas production was close to zero. This 
has been suggested for individual comets - for instance, Levison 
et al. (12006b found that the orbital properties of 2P/Encke may 
be explained, if this comet spent a long time in the past as a 
dormant comet and was relatively recently woken up, as its per¬ 
ihelion distance was decreased by a secular resonance. 

In case dormancy can become permanent, i.e., comets stay 
dormant until they are dynamically ejected, this may influence 
our estimate of past erosion for the reason that the current ac¬ 
tivity of 67P would indicate the comet to be relatively young 
in the dynamical sense. This youth is of course not verified by 
long-term orbital integrations due to the strongly chaotic behav¬ 
ior, i.e., short Lyapunov time scales, of Jupiter family comets. 
Permanent dormancy would seem to be supported by the work 
of Duncan and Levison (11997b . who concluded that the physical 
lifetime of these comets must be limited because their orbits re¬ 
tain a memory of the very low inclinations in the scattered disk 
that they assumed as the source for captures. However, the scat- 
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tered disk is now recognized to be much more excited, so this 
argument may not be of relevance any more. 

We tentatively conclude that a thick surface layer has likely 
been eroded from the 67P nucleus by previous gas-producing ac¬ 
tivity, although the above-mentioned estimate of several hundred 
meters should perhaps be regarded as an upper limit. In any case, 
this cautions that activity-driven erosion may result in complex 
landscapes (Malin and Zimbelman ll9861 Thomas et al. 120051) . 
which could affect our ability to interpret geomorphological fea¬ 
tures in terms of primordial processes. This is particularly true 
for equatorial regions and the southern hemisphere due to their 
computed higher erosion rates (Keller et al. 120151) . Furthermore, 
it should be noted that the small-scale details seen on the sur¬ 
face of the 67P nucleus are likely of very recent origin in most 
cases. For instance, because 67P currently crosses the main as¬ 
teroid belt, following the methodology developed by March! et 
al. (12010112012h .b') we estimate that some 100 - 1000 craters 
larger than a meter should have formed over a time scale of 10^ y 
- however, these craters are not observed. 

For these reasons, we will focus our attention on the over¬ 
all bi-lobate shape and the large-scale structural properties of 
two key regions, Seth and Hathor (Thomas et al. 120151) . which 
are least affected by activity-driven erosion (Keller et al. 2015). 
These regions face the “neck” from the larger and smaller lobes, 
respectively. 

2.2. Contact binary 

The bi-lobate shape of 67P has important implications for the 
formation of the nucleus. Computations based on the current 
shape, spin axis orientation, and a homogeneous composition 
show that the neck-facing regions (in particular, the Seth and 
Hathor regions) are not subject to an increased dose of solar ra¬ 
diation that would imply preferential erosion (Sierks et al. 120151 
Keller et al. l2015t . Although it cannot be ruled out that 67P had 
a radically different spin axis orientation in the past, it seems 
unlikely that erosion of a homogeneous object alone could ex¬ 
plain the current shape of 67P. It is nevertheless possible that 
a very heterogeneous nucleus (with locally higher concentra¬ 
tion of volatiles, or different terrain properties) could result in 
preferential, local erosion. Such a heterogeneous structure, how¬ 
ever, may require mixing of planetesimals with different volatile 
compositions or abundances, probably only attainable through 
collision and merging of objects formed at different heliocentric 
distances. 

Alternatively, the bi-lobate shape may strongly suggest 67P 
to be a contact binary. While contact binaries have been advo¬ 
cated for elongated comet nuclei (e.g., IP/Halley, 19P/Borrelly 
and 103P/Hartley 2), the only unquestionable case remains that 
of 8P/Tuttle based on radar imaging (Harmon at al. 120101) and 
indirectly conhrmed by HST observations (Lamy et al. 120081 ). 
Indeed, based on its appearance, 67P is a highly likely candi¬ 
date. As such, however, its shape appears unusual, since the two 
lobes have their longest axes nearly perpendicular whereas they 
are aligned in most if not all other contact binaries. 

Binary systems are ubiquitous in various small-body popu¬ 
lations, including the near-Earth asteroids (NBAs), main-belt as¬ 
teroids, Jupiter Trojans, Centaurs, and Trans-Neptunian Objects 
(TNOs). Contact binaries alone are estimated to constitute 10 - 
20% of all NBAs, Trojans, and TNOs (Sheppard and Jewitt 
[20041 Noll et al. [20081 Mann et al. [20071 Benner et al. [20061 
12008b . Several mechanisms for the formation of binary systems 
have been proposed, ranging from fission of small rapidly ro¬ 
tating NBAs (Scheeres 12007b . to the formation of large binary 


TNOs by gravitational collapse (Nesvorny et al. 12010b . Whether 
or not these processes are applicable to highly porous, small 
comets deserves further investigation. 

2.3. Structural features 

Concerning possible structural features, the Seth region displays 
multiple semi-planar terrains, or facets (Thomas et al. 120151 
Massironi et al. 12015b . Planar terrains are, in some cases, remi¬ 
niscent of what has been seen on other comet nuclei (e.g., Belton 
et al. 12007b . A key difference, however, is that Seth is charac¬ 
terized by the presence of deep collapsed pits (Vincent et al. 
12015b . whose dust-free walls allow us to glance at the local ver¬ 
tical stratigraphy. In many instances, these walls reveal quasi¬ 
parallel lineaments and small terraces, which may be interpreted 
as the surface expression of internal layering (March! et al. 120151 
Massironi et al. 12015b . 

In the case of the largest terrace in the Seth region, a shape 
model (Jorda et al. 12015b was used to fit the floor of the pit with 
a plane to obtain an approximate local orientation of the puta¬ 
tive layer (Pig.|2]r). Interestingly, the orientation of the resulting 
plane, extrapolated to the extreme end of the large lobe, matches 
the orientation of a topographic ridge used to mark the boundary 
of two very morphologically distinct terrains (Ash and Imhotep 
regions; Thomas et al. 12015b . suggesting that the flat floors and 
lineaments are truly indicative of internal layering, with an esti¬ 
mated local depth of at least a few hundred meters (see March! 
et al. 120151 Massironi et al. 12015b . 

The Seth region faces a prominent dust-free cliff in the small 
lobe, named Hathor, which reveals the ubiquitous presence of 
fractures and terraces, the latter being interpreted as the result 
of layering (Pig. |3^). The two sets of layering in the Seth and 
Hathor regions have local spatial orientations that are not com¬ 
patible with being the expression of a homogeneous layering 
wrapping the two lobes (March! et al. 120151 Massironi et al. 
12015b . A similar conclusion is also reached by looking at the 
global distribution of semi-planar terrains on the surface imaged 
so far (Massironi et al. 12015b . 

The Hathor region also hosts another large-scale geological 
feature with possible implications for the formation of 67P, i.e., 
an impressive set of fractures that run for hundreds of meters 
from the base of the small lobe to its summit (Pig. [3^). The 
largest of these are inferred to be several meters wide, based on 
the cast shadows. Several hypotheses for their formation have 
been investigated, such as thermal fracturing, and gradual sand¬ 
blasting erosion due to activity from the neck region. Thermal 
fracturing typically results in polygonal patterns, such as those 
observed in permafrost on Mars and Earth (e.g.. Bevy et al. 
12010b . and is therefore an unlikely explanation for this well- 
organized set of large fractures. Sand-blasting from a cometary 
jet generally exerts a pressure much lower than the relevant ma¬ 
terial strength, and therefore also does not appear to be a likely 
mechanism. 

Alternatively, wide spread fracturing could be achieved 
through large-scale torques due to re-orientation of the two 
lobes, or by rotational stresses exceeding the failure strength. 
However, such torque-induced fractures are expected to form in 
the neck region itself, where most of the resulting stresses would 
concentrate, rather than the neck-facing regions. Moreover, 
given the almost symmetric nature of the Seth and Hathor re¬ 
gions with respect to the neck, such processes should result in 
similar fracturing. While some degree of fracturing is indeed 
observed in the Seth region, the lack of an extensive and well- 
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organized system of fractures suggests that the two lobes may 
have distinct properties. 

In addition, a spatial analysis of the orientations of the frac¬ 
tures suggests they penetrate deeply into the small lobe, indicat¬ 
ing that they may be pervasive. Support to this conclusion comes 
from the morphology of a large depression, ~1 km across, at the 
summit of the small lobe, named Hatmehit region (Thomas et al. 
12015b . Whether this depression is the result of a collapse or an 
outburst, its rhombic shape indicates the presence of planes of 
weakness within the small lobe, which appear to be aligned with 
the Hathor fracture system (Fig. |3};-d). It therefore seems pos¬ 
sible that these fractures are structural features, pervading the 
smaller lobe. 


2.4. Binary formation 

The above geomorphological evidence on 67P supports the view 
that the bi-lobate nature of 67P is due to collision and merging 
of two distinct objects. Further constraints on this event can be 
derived from the high porosity of 67P (Sierks et al. 12015b and the 
content of volatile species. The high porosity requires peak pres¬ 
sures likely not exceeding ~ lOOkPa (Yasui and Arakawa l2009b . 
and hence an accretional speed not exceeding a few tens of m/s. 
The average impact speed of small planetesimals is thought to 
have evolved from ~ 10 m/s during the early stages still domi¬ 
nated by gas (Kokubo and Ida l2000b to ~ 500 m/s corresponding 
to viscous stirring of the planetesimal disk by its largest mem¬ 
bers (Levison et al. 1201 lb . These considerations narrow down 
two possible formation scenarios for 67P. 

In one case, 67P is the result of a low-speed collision be¬ 
tween two planetesimals, which most likely took place in the ear¬ 
liest stages of the disk evolution before gas dispersal. This would 
mean that the structure of 67P dates back to the time of planetesi¬ 
mal formation and has not been significantly altered by later evo¬ 
lution. Alternatively, 67P is the result of re-accumulation of frag¬ 
ments ejected during a catastrophic collision involving a larger 
planetesimal. Numerical simulations (e.g., Jutzi et al. 12010b have 
shown that the catastrophic disruption of porous parent bodies 
may result in ejection velocities < 10 m/s (comparable to the 
mutual impact speed among fragments), while for a typical im¬ 
pact speed of ~ 500 m/s or lower a significant volume of the par¬ 
ent body may experience low shock pressures (e.g., Leinhardt 
and Stewart [20091 rescaled to lower impact speeds). In addition, 
a highly porous parent body will greatly help dumping the shock 
wave, therefore reducing the pressure exerted far from the im¬ 
pact location (Wiiennemann et al. 12006b . Therefore, the colli- 
sional scenario should also be kept in consideration as a possible 
mode of origin of 67P. In this case too, the structure formed dur¬ 
ing the re-accumulation must not have been destroyed by sub¬ 
sequent evolution, but preceding collisions on the same parent 
body cannot be excluded. 

Our discussion suggests that the two lobes of 67P are two 
distinct objects that merged in a low-speed collision. However, 
the issue remains, whether this collision was part of the primor¬ 
dial planetesimal accretion, or if it involved fragments produced 
later by a collision experienced by a larger parent body. This is 
very important, not only to realize which scenario would best 
explain the properties of 67P as observed by OSIRIS and other 
instruments, but also because current concepts on solar system 
evolution have important implications for the collisional evolu¬ 
tion of comets. We shall now turn to the latter. 


3. Collisional evolution in the trans-planetary disk 

3.1. Generai framework 

We place our discussion within the current paradigm of solar 
system evolution as defined by the Nice Model (Tsiganis et 
al. I20051 for the latest development see Levison et al. I201 lb . 
This model features a disk of planetesimals beyond the initial, 
planetary orbits, which should have extended roughly from 15 
to 30 AU in heliocentric distance. The disk was formed in the 
nascent solar system, when the planetesimals were formed, and 
it was dispersed as a consequence of planet migration following 
the dynamical instability of the giant planets. 

The dispersal is thought to have taken place about 4.1 - 
4.2 Gy ago (Morbidelli et al. 120121 March! et al. 12013b . implying 
~ 0.4 Gy of prior collisional evolution. The total mass of the disk 
is constrained to be in the range 20-50 Earth masses (Tsiganis et 
al. 12005b . The timing of the dynamical instability is assumed to 
coincide with the start of the Late Heavy Bombardment (LHB). 
Among the important effects of the instability, there is also the 
gravitational scattering of objects into the Oort cloud and the 
scattered disk (Brasser and Morbidelli 12013b . We consider that, 
at the typical velocities of collisions in the trans-planetary disk, 
projectiles with diameter (D) larger than 1 km would shatter tar¬ 
gets of the same size as 67P (D ~ 4 km, assuming a past ero¬ 
sion of a 300 m thick layer). This estimate is likely conservative, 
judging from the expected material properties (e.g., Benz and 

AsphaugHUlB 

Concerning the availability of such projectiles, a recent study 
by Johansen et al. (12015b of planetesimal formation in a pebble 
cloud found the cumulative size frequency distribution (SFD) 
to be quite shallow with a power law index a ^ -1.8 for 
20 < D < 200 km, although the applicability of this shallow 
slope down to km-sized objects needs to be further investigated. 
Additional constraints on this size range come from the observed 
trans-Neptunian populations and Trojan asteroids (the current 
leftovers of the primordial disk). These populations exhibit two- 
sloped cumulative SFDs, with a break at D ~ 100 km. The slope 
of the cumulative SFDs for smaller objects ranges from -2.0 to 
-2.5 (Grav et al. l201 11 Fraser et al. 12014b . Of course, if the plan¬ 
etesimal population evolves by collisions, new small objects are 
likely produced as fragments of larger ones. This is guaranteed 
by the above-mentioned shallow slopes of the planetesimal SFD, 
while it would not hold true for a much steeper planetesimal SFD 
of the kind suggested by Belton (12015b . 


3.2. Number of coiiisions in the trans-pianetary disk 

The number of disruptive collisions (A^coii) for 67P can be esti¬ 
mated using a particle-in-a-box formula (Wetherill ll967b . 


A^coll 


NpuTAp 

V 


( 1 ) 


where Np is the number of disk objects, Ap = 7r(Rc + Rp)^ is 
the cross section for a collision between 67P with average radius 
Rc - 2 km and a partner with radius Rp, u is the average rela¬ 
tive velocity, T is the time interval, and V is the volume of the 


^ Morbidelli and Rickman (12015b used the Benz and Asphaug H999b 
scaling law for hard ice, hit at 1 km/s, and found that the minimum di¬ 
ameter of a disruptive projectile for a target of the same size as 67P was 
always about 1 km or less throughout the primordial disk. However, for 
porous targets the minimum projectile diameter would be smaller, using 
the weak ice scaling law of Leinhardt and Stewart ( 12009b . 
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disk. Assuming viscous stiiring of the disk by its largest mem¬ 
bers (Levison et al. 1201 H . we estimate u - 0.6 km/s, which is 
10% of the typical orbital velocity in the disk. The volume is 
obtained by considering a circular annulus with inner and outer 
radii of 15 and 30 AU, and a constant thickness of 4 AU. 

As collision partners we consider disk objects with radius 
from 0.5 km up to the break point in the cumulative planetesimal 
SFD, which occurs at Rp ~ 50 km. We take a power law with 
constant index a for the cumulative planetesimal SFD and derive 
the corresponding differential number n{Rp)dRp of disk objects 
in a radius range of width dRp. We replace NpAp in Eq. (1) by 
the integral 

^{.^p)ApdRp , 

taking the limits at 0.5 km and 50 km. The computation of n{Rp) 
will hrst be done as follows. At the end of the time inter¬ 
val {T - 0.4 Gy) we use an estimate of 2 ■ 10" objects with 
D > 2.3 km in order to populate the scattered disk (Brasser 
and Morbidelli I2013t . Populating the Oort cloud would seem 
to require even more, but we prefer to be conservative. For the 
cumulative SFD index a we consider two values: -2.5, which 
would correspond to a relaxed, collisional asteroid population 
(Dohnanvi 1 19691 1, and -2.0, which is close to the value found 
by a theoretical model (Johansen et al. 12015t for small planetes- 
imals formed by pebble accretion. The latter would likely char¬ 
acterize the initial state of the disk, and the former may charac¬ 
terize the hnal state after collisional evolution. 

For the sake of illustration, we use the above value for the 
cumulative number with D > 2.3 km with the two a values to 
characterize the disk throughout the time interval, even though 
it explicitly refers to the hnal state. Our calculation then yields 
A^coii - 11 and 9 for a - -2.5 and -2.0, respectively. In that 
case, 67P-sized objects would undergo about ten collisions on 
the average, and the chance to survive without any collision, as 
estimated by a Poisson distribution, is less than 10“^. Thus, for 
67P to survive undisrupted, the initial Np would have to be at 
least a factor 10 smaller than we have assumed. Hence, the hnal 
objects would have to be mostly collisional fragments, and it is 
statistically unlikely for 67P to be among the initial survivors. 

Interestingly, the disk mass in the considered size range, as¬ 
suming a density of 1000 kg/m^ (Brown l20131 l. comes out as 6 or 
18.5 Earth masses {Me) for the two a values. We now consider 
the contribution by objects with D > 100 km. The slope of the 
cumulative SEDs for the trans-Neptunianpopulations and Trojan 
asteroids for diameters larger than ~ 100km ranges from -3.5 
to -4.5. We will use a - —A for the larger objects (Gladman et 
al. 1200 11 Eraser et al. l2014l l while keeping both limits discussed 
above for the smaller size range as in the previous computation. 
Very large objects (such as Pluto) are insufficiently sampled by 
the adopted planetesimal SED, but it is expected that ~ 1000 
Pluto-sized objects once existed in the primordial disk (Levison 
et al. l201 II) . thus contributing a few Earth masses. With these as¬ 
sumptions, the total mass is compatible with the inferred range 
from the Nice Model as stated above. 

An alternative way of estimating n{Rp) is as follows. We 
take the minimum value of 20 Me for the total disk mass. We 
subtract 2 Me for the Pluto-sized objects and consider 18 Me to 
be partitioned between two populations with D < 100 km and 
D > 100km. We distribute the IS Me up to D = 1000km, and 
as a result we find for a = -2.5 that the number of objects with 
D > 2.3 km is 3 ■ lO", yielding Vcoii - 18, and for a - -2.0 
we hnd 7 • 10^°, yielding Vcoii - 3. In this case, even for the 
shallower cumulative planetesimal SED, the chance for survival 


without any collision is only 5%, and yet the planetesimal popu¬ 
lation has an insufficient number of objects to feed the scattered 
disk. Figure|4]illustrates the mass distribution used for the trans¬ 
planetary disk objects in this analysis. 


4. Escape from collisions 

We conclude from the above that in both scenarios under con¬ 
sideration, 67P is likely to be collisionally evolved. To judge the 
chances for 67P to be unaffected and primordial, we may con¬ 
sider a possible gradient of collision probability across the disk. 
Its surface density is usually modeled to vary as with he¬ 
liocentric distance r. In addition, the thickness of the disk may 
have been proportional to r, so that the number density may have 
varied as This would cause Vcoii to be four times lower 
in the outermost parts of the disk than in the innermost parts. 
Therefore, 67P may possibly be a primordial survivor without 
being too exceptional, if it formed in the outskirts of the disk, as 
possibly indicated by its high D/H ratio (Altwegg et al. I2015t . 
In addition, an even lower mass of the disk (Nesvorny, pers. 
comm.), and a less dynamically excited disk could also reduce 
the estimated collision rate by a few times. In any case, it appears 
that the best chance for 67P to be primordial is with a low-mass 
disk and a shallow size distribution of the planetesimals (namely, 
O' < -2 for D < 100km). 

Our model features a common value of the encounter veloc¬ 
ity for the trans-planetary disk objects. This is thought to be an 
average of the individual velocities, relevant to the average state 
of dynamical excitation of the disk. Of course, there are many 
encounters occurring at both lower and higher velocities, so we 
should consider the effect of this scatter. Since we used a cut-off 
of projectile diameters at D = 1 km, we cannot account for the 
fact that the size limit for disruptive collisions decreases with in¬ 
creasing velocity. However, this together with the fact that the 
collision frequency increases with velocity means that we have 
likely underestimated the number of disruptive collisions. 

Another caveat concerns the timing of the giant planet in¬ 
stability. Even though it is considered likely that the LHB was 
triggered by the flux of asteroids and comets destabilized by the 
planet migration, there remains a possibility that the instability 
occurred earlier than assumed here. In such a case, the disk life¬ 
time may have been signihcantly less than 400 My, thus dras¬ 
tically lowering the above estimates of Ncoii. However, a very 
early instability would face other problems such as hnding an 
alternate explanation for the LHB. 

We should also keep in mind that, whatever was the fate of 
the planetesimals from which 67P stems, our estimates confirm 
the recent conclusion by Morbidelli and Rickman (12015b that 
comets in general did not escape collisional evolution in the early 
solar system. This holds for all kinds of comets - those coming 
from the Oort cloud as well as those captured from the scattered 
disk or the EKB. Therefore, if 67P and other comets are found to 
retain primordial structural features, it follows that these cannot 
have been erased by the collisional history of the comets. 

In this perspective, 67P appears to offer a precious test case, 
since it has not stood out as very special among the JECs. Thus, 
the study of a collisional re-accretion origin for 67P is partic¬ 
ularly important. For instance, should such an origin prove to 
be inconsistent with the body of evidence from OSIRIS and 
other Rosetta instruments in the sense that the primordial sig¬ 
natures {e.g., layering and outgassing of super-volatiles) would 
have been erased, this would suggest that both 67P and comets at 
large, are indeed primordial survivors. Such a conclusion would 
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necessitate a reconsideration of the concepts around the early 
evolution of the solar system that were applied in this paper. 

On the other hand, should it turn out that collisional re¬ 
accretion is consistent with the observations of 67P including 
those features thought to pre-date this accretion (e.g., layering 
and volatile abundances), the study of 67P would still reveal im¬ 
portant clues about early planetesimal formation, independent of 
the role played by collisions in its later history. 


5. Discussion and conciusions 

We have shown that OSIRIS data helps to constrain the origin of 
comet 67P/Churyumov-Gerasimenko in several important ways. 
The bi-lobate overall shape is best interpreted in terms of a con¬ 
tact binary, formed by the collisional merger of two distinct bod¬ 
ies. Their properties appear generally similar, but the smaller 
lobe may possibly show evidence of damage due to an energetic 
impact, perhaps witnessed by the observed large-scale fractures. 

The large porosity derived from the bulk density of 
470 kg/m^ in fact requires accretional velocities for the two lobes 
< 10 - 20 m/s, and this leaves open two options for the origin 
of 67P. Either it may be an undestroyed, primordial planetesi¬ 
mal that accreted during the hrst ~ 10 My, or it was formed by 
the gentle re-accretion of fragments ejected from a larger par¬ 
ent body upon a destructive collision at a later time (possibly, 
the last one of a series of collisions). It is not possible yet to ex¬ 
clude either of these scenarios. However, should the latter option 
be true, it must have conserved the primordial signatures of 67P 
(both geomorphological and chemical), and hence the comet re¬ 
mains a key witness of planetesimal formation in the nascent 
solar system. We also emphasize that the collisional scenario re¬ 
quires the bi-lobate nucleus to have survived intact after it was 
formed, and thus the collision that led to its formation must have 
occurred during the hnal stages of the primordial disk or even 
during its dispersal. 

The origin of internal layers remains elusive. However, it 
should be noted that layering appears to be an inherent char¬ 
acteristic of the two lobes prior to their accretion. According to 
Groussin et al. (120151 1. the low compressional strength of 67P 
may stimulate pressure-induced layering. The early, gentle dock¬ 
ing of two km-sized planetesimals thus implies that their layer¬ 
ing arose during their primordial accretion. Alternatively, the in¬ 
volvement of a catastrophic disruption of a larger parent body (a 
few 10s km at most to avoid signihcant reduction of porosity by 
compaction) could explain layering, if some internal evolution 
had taken place resulting in a onion-shell structure, and the two 
lobes are large chunks of the parent body rather than debris piles. 

We also show that strong collisional evolution is generally 
to be expected for the objects that occupied the early, trans¬ 
planetary disk. While it is possible for 67P to have survived 
collisional destruction in a safe niche under special assumptions 
about the disk, we point out that 67P has not shown any evidence 
of being exceptional among comets at large, and thus it should 
not have an exceptional origin. However, this also means a chal¬ 
lenge for the collisional origin scenario. Low densities have been 
found for all Jupiter family comets so far analyzed, and if these 
are generally produced by destructive collisions involving large 
parent bodies, the latter must have conserved deep surface layers 
in a primordial, porous state in spite of Al-induced chemical 
differentiation (Prialnik et al. I2004I I. Interestingly, several EKB 
objects in the size range D ~ 150 - 400 km have been found 
to have a density of 400 - 800kg/m^ (Brown 12013l l. implying 
signihcant porosity can be retained even on large objects. 


Regarding the two main options, we conclude that the pri¬ 
mordial survivor option would likely call into question some of 
the current concepts of solar system evolution, while the colli¬ 
sional fragment option would not call for any such rethinking. 
The issue with the latter is to show it to be consistent with the 
observational evidence on 67P and other comets. 
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Fig. 1. Evolutions of the perihelion distance of 67P, found by forward integrations of random initial orbits leading into current orbits 
close to the one of 67P. The time plotted on the abscissa is counted backward from the present. The different colors mark the different 
examples chosen. We note that these are not clones of 67P, and the recent history of 67P (including the close encounter with Jupiter 
in 1959) is not reproduced. The focus is on the long-term dynamics, shown in the right part of the diagram. Chaotic divergence of 
the backward motions sets in, mostly a few centuries ago, due to close encounters with Jupiter. Courtesy Ryszard Gabryszewski. 
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Fig. 2. Large-scale layering on the large lobe, a) Close-in view of the pitted terrain on the large lobe (Seth region). The red lines 
indicate some of the prominent lineaments, interpreted as surface expression of layering, b) A plane indicating the prevailing local 
orientation of the layers for the largest pit in the Seth region, c) A set of lineaments and high topographic reliefs on the opposite side 
of the large lobe, d) Same plane as in b) and two additional parallel planes. The orientation of these planes matches the lineaments 
and topographic reliefs, suggesting the latter are surface expressions of internal layering. We note that the local orientation of the 
plane may vary from place to place. For simplicity we assumed that the best fit plane does not change attitude within the large lobe. 
However, a similar conclusion applies for bending following the local convex shape, or along the y-axis in panel b. 
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Fig. 3. Large-scale fractures on the small lobe, a) Cross-cutting pattern of fractures (orange lines) and layers (purple lines) on the 
Hathor region, b) A plane indicating the prevailing orientation of the fractures overlaid on a shape model, c) A view of the top of the 
small lobe, with some major features highlighted: Layering on the large lobe (red lines; see Fig. 13, rim of the Hatmehit depression 
(yellow lines), d) The same plane shown in b), and an additional parallel plane. The orientation of these planes matches the direction 
of two sides of the Hatmehit depression. 
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Fig. 4. A sketch of the primordial disk cumulative mass distribution (in units of Earth’s masses), based on dynamical constraints 
and current observational constraints from Kuiper Belt objects and Trojan asteroids. We plot, as a function of diameter, the total 
disk mass contained in all larger objects. Two Earth masses have been allocated to objects larger than 1000km. The total mass of 
objects larger than 1 km has been set to 20 Earth masses, at the low end of published primordial disk masses. The two solid curves 
correspond to slopes of -2.0 and -2.5 for the cumulative size-frequency distribution at D < 100 km. 
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